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ABSTRACT 

We present a if^-selected catalog, dubbed FIREWORKS, for the Chandra Deep Field South (CDFS) 
containing photometry in C/38, B435, B, V, Veoe, R, «775, I, ^850, J, H, Kg, [3.6/im], [4.5/im], [5.8/im], 
[8.0/im], and the MIPS [24/iin] band. The imaging has a typical i^*°AB limit of 24.3 mag (5cr) and cov- 
erage over 113 arcmin^ in all bands and 138 arcmin^ in all bands but H . We cross-correlate our catalog 
with the 1 Ms X-ray catalog by Giacconi et al. (2002) and with all available spectroscopic redshifts 
to date. We find and explain systematic differences in a comparison with the '^sso + if s '-selected 
GOODS-MUSIC catalog that covers ~90% of the field. We exploit the C/38-to-24 ^m photometry to 
determine which i^s-selected galaxies at 1.5 < z < 2.5 have the brightest total IR luminosities and 
which galaxies contribute most to the integrated total IR emission. The answer to both questions 
is that red galaxies are dominating in the IR. This is true no matter whether color is defined in the 
rest-frame UV, optical, or optical-to-NIR. We do find however that among the reddest galaxies in 
the rest-frame optical, there is a population of sources with only little mid-IR emission, suggesting a 
quiescent nature. 
Subject headings: galaxies: distances and redshifts - galaxies: high rcdshift - infrared: galaxies 



1. INTRODUCTION 

Since the original Hubble Deep Field (Williams et 
al. 1996), deep multi-wavelength observations of blank 
fields have revolutionized our understanding of the high- 
rcdshift universe. Especially the epoch around z ~ 2 
is of great interest since it is then that the cosmic star 
formation rate density was peaking (Hopkins & Beacom 
2006). At z ^ 2, the observed optical probes the red- 
shifted rest-frame UV emission, making it a good tracer 
for relatively unobscured star formation. Near-Infrared 
(NIR) observations of distant galaxies, such as under- 
taken by the FIRES survey in the Hubble Deep Field 
South (HDFS, Labbe et al. 2003, hereafter L03) and the 
MS 1054-03 field (Forster Schreiber et al. 2006, hereafter 
FS06), show relatively small variations in the rest-frame 
V^-band mass-to-light ratio. Selecting galaxies in the Kg- 
band (e.g., L03; FS06) thus provides a good probe of the 
massive galaxy content at high redshift. 

In the presence of dust, large amounts of rest-frame 
UV emission can be absorbed and re-emitted in the Far- 
Infrared (FIR) . Dust corrections of the UV luminosities 
of such systems involve large uncertainties. Direct obser- 
vations of the dust emission are therefore crucial to get a 
better census of the bolometric energy output. In order 
to study the bolometric properties of typical galaxies at 
2 ~ 2, infrared luminosities have been derived from the 
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observed 24 /im flux by means of IR spectral energy dis- 
tribution (SED) templates (e.g., Papovich et al. 2005; 
Reddy et al. 2006). 

Despite the extra model uncertainty involved, this ap- 
proach adds complementary information to the shorter 
wavelength studies of high-redshift galaxies. While the 
Luminous (lO^^L© < Lir < lO^'^^i©) and Ultraluminous 
(LiR > IQI^Lq) Infrared Galaxies, (U)LIRGS (Sanders 
& Mirabel 1996) are locally very rare, they are found 
to be increasingly more common toward higher redshifts 
(e.g., Caputi et al. 2006). 

In this paper, we present a Kg-hand selected multi- 
wavelength catalog for the GOODS-CDFS, comprising 

WFI UssBVRI, ACS B435^606*7752850, ISAAC JHKs, 

IRAC 3.6-8.0 ^m and MIPS 24 ^m imaging. We adopt 
a similar format as for the publicly available FIRES cat- 
alogs of the HDFS and MS 1054-03, hence the name 
FIREWORKS. This allows the user to exploit the com- 
bined photometry of the CDFS, MS 1054-03, and the 
HDFS in a straightforward manner. The fields are com- 
plementary in depth (5(T for point sources Kl°^^ = 24.3, 
KI°Xj^ = 25.0, X*°|b = 25.6 respetively) and area (138 
arcmin^, 24 arcmin^, and 5 arcmin^ respectively). 

An analysis of the space density and colors of massive 
galaxies at 2 < z < 3 (van Dokkum et al. 2006) , the rest- 
frame optical luminosity density and stellar mass density 
up to z ~ 3 (Rudnick et al. 2006), and the rest- frame 
luminosity functions of galaxies at 2 < z < 3.5 (March- 
esini et al. 2006) were partly based on the FIREWORKS 
catalog for the GOODS-CDFS presented here. 

After describing the catalog construction, we partic- 
ularly address the questions which A's-selected galax- 
ies at 1.5 < z < 2.5 are brightest at 24 fim, which 
galaxies have the largest total infrared luminosity Lir, 
(= L(8 — 1000 /im)) and contribute most to the total in- 
tegrated IR luminosity emitted by Kg-selected galaxies. 



We address this question by studying the IR emission 
as function of color defined in three wavelength regimes: 
the rest-frame UV, optical, and optical-to-NIR. 

An overview of the observations is presented in §2. §3 
describes the construction of the final mosaics. Source 
detection and photometry is discussed in §4. Next, 
we present our photometric redshifts (zphot) and cross- 
correlation with the available spectroscopic surveys in 
§5. §6 summarizes the catalog content. A photometric 
comparison for the wavelength bands in common with 
the GOODS-MUSIC catalog by Grazian et al. (2006a, 
hereafter G06a) and a ^phot comparison with the same 
authors is discussed in §7. Results on 24 ^m properties 
and total infrared luminosities of iiTs-selected galaxies at 
1.5 < z < 2.5 are discussed in §8. §9 summarizes the 
paper. 

AB magnitudes are used throughout this paper. We 
adopt a iJo = 70 km s'^ Mpc~i, flM = 0.3, 17a = 0.7 
cosmology. 

2. OBSERVATIONS 

2.1. The GOODS Chandra Deep Field South 

Centered on {a,S) = (03:32:30, -27:48:30), the CDFS 
(Giaconni et al. 2000) has been targeted by most of 
today's major telescope facilities, both in imaging mode 
over the whole spectral range and in spectroscopic mode. 
Among the surveys carried out in areas included in the 
CDFS are K20 (Cimatti et al. 2002), GMASS (Cimatti 
et al. in preparation), VVDS (Le Fevre et al. 2004), 
and the Hubble Ultra Deep Field (Beckwith et al. 2006). 
In this section, we describe the public GOODS-South 
datasct that we used to build a Xg-band selected catalog 
containing homogeneous colors from the optical to 24 
/xm. 

2.2. The WFI UasBVRI data 

Ground-based optical imaging of the Extended CDFS 
with the Wide Field Imager (WFI) on the ESO/MPG 
2.2m telescope was obtained as part of the COMBO-17 
(Wolf et al. 2004) and ESO DPS (Arnouts et al. 2001) 
surveys. We use the central part of the WFI U^s, B, V, 
R, and / imaging, where it overlaps with the GOODS 
Ks-hand imaging, that has been re-reduced by the Ga- 
BoDS consortium (Erben et al. 2005; Hildebrandt et al. 
2006). The images were released on a 0'.'238 pixel^^ scale 
and have exposure times of 49, 69, 105, 88, and 35 ks in 
C/38, B, V, R, and / respectively. 

2.3. The ACS B435Veoei77zZs50 data 

During 5 epochs of observations, the ACS camera on 
HST acquired imaging of the GOODS-South field in 4 fil- 
ter bands: F435W, F606W, F775W, and F850LP (here- 
after referred to as i?435,V606:*775j and zsbo)- Exposure 
times amounted to 2, 1.7, 1.7, and 3.3 hours respectively. 
The mosaics (version vl.O; Giavalisco et al. 2004), were 
drizzled onto a pixel scale of 0'.'03 pixel" ^. From the 
150 arcmin^ area that is well covered by the Ks-ha,nd 
detection image, 138 arcmin^ is well exposed with ACS. 
We restrict our analysis of the total IR properties of the 
distant galaxy population in §8 to this overlap region. 



We use the ESO/GOODS data release vl.5 
(http://www.eso.org/science/goods/releases/20050930/) 
to complement the optical observations with NIR imag- 
ing by the Very Large Telescope (VLT). For a full 
description of the dataset, we refer to Vandame et al. 
(in preparation). Briefly, the vl.5 data release consists 
of 24 fully reduced VLT/ISAAC fields in the J and K^ 
bands and 19 fields in the -ff-band, each with a 2f5 x 2.'5 
FOV and 0'.'15 pixel" ^ scale. The ISAAC data were 
reduced using the ESO/MVM image processing pipeline 
(vl.9, see Vandame 2002 for the description of an earlier 
version). Exposure times varied from field to field, with 
typical exposures of 3.2, 4.2, and 5 hours in J, H, and 
iiTs-band respectively, and respective ranges between 
ISAAC fields of 2-5 hours, 2-6 hours, and 3.6-7.5 hours. 
The variations in depth resulting from the unequal 
exposure times are discussed in §3.5. A total area of 
113 arcmin^ is well exposed in all optical and NIR 
filter bands. Without a restriction on the iJ-band, the 
covered area increases to 138 arcmin^. 

2.5. The IRAC 3.6-8.0 fim data 

As a Spitzer Space Telescope Legacy Program, su- 
perdeep images of the GOODS-South field were taken 
with the Infrared Array Camera (IRAC, Fazio et al. 
2004) on-board Spitzer. Over 2 epochs the whole field 
was covered in the 3.6 /xm, 4.5 fim, 5.8 /xm, and 8.0 /ini 
bands. For each epoch, exposure times per channel per 
sky pointing amounted to 23 hours. With the telescope 
orientation being rotated by 180 degrees between the two 
epochs, the second epoch IRAC channel 1 and 3 observa- 
tions targeted the area covered by IRAC channel 2 and 4 
during the first epoch, and vice versa. An overlap region 
of roughly 40 arcmin^, including the Hubble Ultra Deep 
Field (Beckwith et al. 2003), got twice the exposure time. 
We use the data releases DR2 and DR3 for the second 
and first epoch respectively. Images were released on a 
0'.'60 pixel" '^ scale. A full description of the observations 
and reduction will be presented by Dickinson et al. (in 
preparation). 

2.6. The MIPS 24 fim data 

The GOODS-South field was observed at 24 /im with 
the Multiband Imaging Photometry for Spitzer (MIPS, 
Rieke et al. 2004) on-board Spitzer, closely overlapping 
the IRAC fields with a position angle that is rotated 
with respect to the IRAC observations by approximately 
3 degrees. The MIPS campaign led to a nearly uniform 
exposure time of 10 hours. We use the version v0.30 
reduced images, released on a 1'.'20 pixel" ^ scale, based 
on the Spitzer Science Center (SSC) Basic Calibrated 
Data (BCD) pipeline (version Sll.0.2). 

3. FINAL IMAGES 

In this section, we describe the image quality of the 
publicly released data products, the subsequent steps un- 
dertaken to obtain the final mosaics from which the pho- 
tometric catalog is extracted, and the limiting depths 
reached at all wavelengths. The main characteristics of 
the multi-wavelength data set are summarized in Table 
1. 



2.4. The ISAAC JHK, data 



3.1. Pixel scales and large scale backgrounds 



Table 1. Summary of the multi- wavelength data set for the GOODS- 
south field 

Camera Filter Ac FWHM^ Sct Total Limiting Magnitude 

(A) (arcsec) (AB mag) 



Ground-based Optical Bands 



WFI 
WFI 
WFI 

WFI 
WFI 



B 
V 
R 
I 



3656 
4601 
5379 
6516 

8658 



0.93 
0.97 
0.90 
0.78 
0.93 



26.06 
27.21 
26.88 
26.99 
25.00 



Space-based Optical Bands 



ACS B43B 

ACS Veoe 

ACS J775 

ACS zsso 



4327 
5957 
7705 
9072 



0.22 
0.22 
0.21 
0.22 



27.68 
27.81 
27.26 
26.90 



Near-infrared Bands 



ISAAC J 



12379 0.38 - 0.61 



25.82 - 26.42 



ISAAC 
ISAAC 



H 



16517 0.35 
21681 0.38 



0.65 
0.58 



25.39 ■ 
25.02 ■ 



Mid-infrared Bands 



IRAC 


chl 


35634 


1.6 


IRAC 


ch2 


45110 


1.7 


IRAC 


ch3 


57593 


1.9 


IRAC 


ch4 


79595 


2.0 


MIPS 


24 Atm 


237000 


6.0 



25.95 
25.96 



26.15 
25.66 
23.79 
23.70 
21.30 



''The intrinsic ACS PSF has a FWHM of 0'.'12. The values quoted in this table were measured on the 5x5 block averaged mosaics. The 
ACS photometry was performed on the images convolved to the broadest NIR PSF. 



First, we converted the ACS images to the 0'.'15 pixel"^ 
scale of the Xg-band detection image by 5 x 5 pixel block 
averaging with flux conservation. All space-based optical 
and NIR photometry was performed on this pixel scale 
using the SExtractor software version 2.2.2 (Bertin & 
Arnouts 1996) (see §4). 

A source fitting algorithm developed by Labbe et al. 
(in preparation), especially suited for heavily confused 
images for which a higher resolution prior (in this case 
the Xg-band image) is available, was used to extract pho- 
tometry from the WFI, IRAC and MIPS images (see 
§4.2.0). The WFI images were registered to the i^T^-band 
mosaic with a 0'.'15 pixel" ^ scale using position measure- 
ments of stars and compact objects and the IRAF ge- 
omap and geotran tasks. For the MIR photometry, the 
algorithm requires a higher resolution image than pro- 
vided by the IRAC and MIPS images. However, the na- 
tive iiTs-band pixel scale makes it more computationally 
expensive without benefit in accuracy. A 2 x 2 block av- 
eraged version of the i^Tg-band mosaic was therefore pro- 
duced with a O'.'S pixel" ^ scale. We registered the pub- 
licly released IRAC and MIPS images onto this iiT^-band 
image using the WCS information in the image header in 
combination with a minor additional shift based on po- 
sition measurements of stars and compact sources, again 
forcing fiux conservation. We note that the source fit- 
ting algorithm takes care of residual shifts. Since the 
programme does not take into account large scale back- 
ground variations, these were removed a priori from the 



IRAC and MIPS images by subtracting SExtractor back- 
ground images produced with large background mesh 
settings (BACK_SIZE=64, BACK.FILTERSIZE=3 on a 
O'.'S pixel" ^ scale). The subtracted background varied 
over the image by -2 to times the rms noise. 

3.2. Image quality and PSF matching 
Space-based optical and NIR wavelengths 

In order to obtain consistent color measurements, we 
match all space-based optical and NIR images to a com- 
mon resolution, namely that of the field with the broad- 
est point spread function (PSF). Since the ground-based 
optical, IRAC, and MIPS images have a significantly 
broader PSF, they are not included in the simple PSF 
matching, but are treated separately by a source-fitting 
algorithm (see §4.2.0). In this section, we describe the se- 
lection of stars used to build the PSFs, the construction 
of the PSFs for the ACS -8435-, Veoe-, «775-, and zsso- 
mosaics and for each of the ISAAC fields in the J-, H-, 
and i^s-band, the construction of the convolving kernels, 
and the quality of the PSF matching. 

First, we compiled a list of bright, isolated, unsatu- 
rated stars. Initially, well-exposed objects with (J — 
Ks) AB < 0.04 and if *°ab < 22.86 mag were selected from 
a preliminary catalog of the CDFS. For ISAAC fields 
where the number oi J ~ Kg selected stars was low, we 
complemented the sample with stars from the EIS stellar 
catalog (Grocnewegen et al. 2002). During a first itera- 
tion, the list was cleaned from galaxy-like objects, stars 
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Fig. 1.— Distributions of seeing FWHM for the ISAAC J, 
H, and Ks observations. Moffat profiles were fitted to the 
PSFs that were built from bright, isolated, unsaturated stars 
for each field and band separately. 



with neighbors within 3" radius, stars too close to the 
edge of an image, objects that were not identified in the 
ACS rl.lz catalog (Giavalisco et al. 2004) or resolved 
objects as defined by having a FWHM in the zsso-band 
larger than O'.'IS. Measurements of the FWHM were per- 
formed by fitting Moffat profiles to the stars using the 
imexam task in IRAF. Finally, we inspect by eye the 
radial profiles and curves of growth, produced with the 
IRAF tasks radprof and phot respectively. The final list 
contained 3 to 5 stars per ISAAC field, with the excep- 
tion of field f30, for which only 1 good star was available. 
The numbers of stars used to build the ACS PSFs were 
31, 45, 49 and 53 for the ^435, Vgoej ^775: and 2:850 mo- 
saics respectively. 

It is sufficient to construct one PSF image per band 
for the whole ACS mosaic since for space-based imaging 
the PSF is uniform across all pointings. The ground- 
based data of each ISAAC field were taken under dif- 
ferent average seeing conditions, hence the need for the 
construction of an individual PSF image for each field 
and each NIR band. We computed the J-, H-, and Kg- 
band PSF images per ISAAC field by averaging the reg- 
istered and flux-normalized images of the selected stars. 
The flux was normalized within 1'.'5 diameter apertures 
rather than the total aperture to optimize the signal-to- 
noise and avoid contributions from residual neighboring 
sources. Any neighbors in the ISAAC images of good 
stars, sufficiently far not to bias the FWHMs and PSFs, 
were masked while averaging. This method was pre- 
ferred over taking the median, since only a handful of 
good stars per ISAAC field were available to build a PSF. 
PSFs for the ACS mosaics were constructed from a large 
enough number of stars to average out any influence of 
faint neighboring sources without masking. 

The ACS PSFs in the -B435-, Vgoe-j *775-: and zgso-band 
on the 5x5 block-avcragcd ACS mosaics are well fit by 
Moffat profiles with /3 = 3.8, 4.0, 3.7, and 3.6, and a 
FWHM of 0'.'22, 0'.'22, 0'.'21, and 0'.'22 respectively. The 



FWHM of the NIR ISAAC observations varied from 0'.'35 
to 0'.'65, with median values of 0'.'47, 0'.'48, and 0'.'47 in 
J, H, and Kg respectively. The values of the Moffat f3 
parameter were 3.1, 2.9, and 3.0 in J, H, and Kg. Fig. 
1 illustrates the distribution of FWHMs of the natural 
PSFs for the individual ISAAC fields. In all of the consid- 
ered bands and fields, the FWHM of the individual stars 
were within «10% of that measured on the final average 
PSF of the image. We adopted the 0'.'65 i?-band PSF of 
ISAAC field fl5 as target to which all higher resolution 
images were matched. 

We computed the kernel for convolution for each 
ISAAC field and band separately, using the Lucy- 
Richardson deconvolution algorithm. The kernels are 
well represented by a Moffat profile with /3 parameters 
varying from field to field and band to band. The ratio 
of the growth curve of the convolved PSF over that of 
the target PSF is a good measure for the PSF matching 
accuracy. Overall, the ratio of growth curves deviated 
by at most 5.1% from unity for apertures between 1" 
(«1.5FWHM of the PSF of the smoothed field maps) 
and 6" (the reference aperture for photometric calibra- 
tion), with an average of 0.54% ± 0.90%. Flux within 3" 
radius is well conserved during the convolution process, 
with a median deviation of 0.06% and a median absolute 
deviation of 0.6%. 

The construction of convolving kernels for the ACS 
mosaics required an extra step. The kernels obtained 
from deconvolution with the IRAF Lucy task had sig- 
nificant noise in the outer parts, leading to noise spikes 
around bright point-like sources in the convolved ACS 
mosaics. To remove these artifacts, we modeled the ACS- 
to-ISAAC kernels by fitting isophotes using the IRAF 
tasks ellipse and bmodel, and used the modeled kernels 
for the convolution. This is possible because the kernels 
are otherwise well behaved and very azimutally symmet- 
ric. The ACS-to-ISAAC kernels have a nearly Moffat 
profile with a /3 parameter of 3. 

Because of the different basic shapes of the ACS and 
ISAAC PSFs, an excellent matching over the relevant 
radii is more difficult than in general among ISAAC 
fields. Nevertheless, the offsets of the growth curve ratios 
between 1" and 6" are limited to below 4.7%, with an av- 
erage of 1.58% ± 1.32%. The average over all stars of the 
ratio of the fiux measured within 3" radius in the con- 
volved and that measured in the natural image showed 
a flux conserving accuracy of 0.7% or better for all ACS 
bands. 

Ground-based optical and MIR wavelengths 

The seeing PSF of the WFI observations is signifi- 
cantly broader than that of our Kg-hemd detection image 
(- 0.9" FWHM). The FWHM measured on the average 
image of bright, isolated stars in the IRAC images even 
amounts to 1'.'6, 1'.'7, 1'.'9, and 2'.'0 for the 3.6 ^m, 4.5 ^m, 
5.8 ^m, and 8.0 ^m bands respectively. The MIPS 24 /im 
beam has a FWHM as large as 6" . Since confusion and 
blending effects are unavoidable in deep observations at 
this resolution, we decide not to degrade the space-based 
optical and NIR images to the MIR resolution. Instead, 
we construct PSFs and convolving kernels similarly as de- 
scribed in §3.2.0, but apply them using a source fitting 
algorithm that makes fully use of the higher resolution 
information in the ii'^-band detection image (see §4.2.0, 



Table 2. H-band zero points in the AB system derived from the NIR 
stellar locus 



Field 


//-band zero point 


03 


25.99 


04 


26.02 


05 


26.07 


08 


25.89 


09 


25.92 


10 


25.94 


11 


25.93 


13 


26.02 


14 


25.82 


15 


26.07 


16 


25.97 


19 


25.86 


20 


25.89 


21 


25.97 


22 


26.03 


23 


26.07 


24 


25.95 


25n 


25.94 


26n 


26.09 



Labbe et al. in preparation). 

3.3. Zero points 

The zero-point calibrations for all bands but the H- 
band were taken from the respective GOODS data re- 
lease. In the case of the NIR ISAAC observations, the 
publicly released zero points were based on SOFI images 
of the EIS-DEEP and DPS infrared surveys conducted 
over the same region (Vandame et al. 2001), which 
themselves were photometrically calibrated using stan- 
dard stars from Persson et al. (1998). That procedure 
yielded zero points with rms scatters ranging between 
0.01 and 0.06 mag in the J-band, 0.01 and 0.08 mag in 
the Ks-ha,nd and up to 0.17 mag in the i?-band. 

To improve on the H-h&nd calibration, we make use 
of stellar photometry in the FIRES HDFS (LOS) and 
MS 1054-03 (FS06) fields, for which i?-band zero points 
were determined to a ~ 0.03 mag accuracy. For each of 
the 19 ISAAC fields with H-hand coverage, we measured 
the mean offset of the stars used for PSF matching along 
the J — H axis of a J — Ks versus J — H color-color dia- 
gram with respect to the stellar locus in the FIRES fields. 
Assuming the J- and ii's-bands are well calibrated, this 
immediately provides us with the if-band zero point cor- 
rections to be applied. We list the derived iJ-band zero 
points in Table 2. Zero-point corrections ranged from - 
0.18 mag to 0.09 mag, with a median correction over all 
fields of -0.03 mag. After applying the zero-point cor- 
rection, the median absolute deviation in J — if color of 
individual stars around the stellar locus is 0.03 mag, sim- 
ilar as measured for the FIRES HDFS and MS 1054-03 
fields. 

3.4. Mosaicing and astrometry 

Here, we describe the combination of the smoothed 
ISAAC NIR fields and the astrometric precision of the 
final mosaics. The 5x5 blocked and smoothed ACS 
Z775-band mosaic was adopted as astrometric reference 
image. The astrometric solution for the ACS data it- 
self was based on a cross-identification of sources with 
deep ground-based WFI data that on its turn was astro- 
metrically matched to stellar positions in the Guide Star 
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Fig. 2. — Map of residual shifts of compact sources in the 
T^s-band mosaic with respect to the reference ACS irrs-band 
mosaic. 2(7-clipped reference sources used for the alignment 
are indicated in black. Grey vectors represent the residual 
shift of the 2(7 outliers. 



Catalog 2 (GSC2, STScI 2001). The final solution had a 
2cr-clipped rms deviation of < O'.'Ol in ACS-to-ACS and 
0'.'12 in ACS-to-ground difference. 

The smoothed ISAAC fields were registered onto the 
smoothed ACS Z775-band mosaic by applying simple x- 
and y- shifts without further distortion corrections. In 
each ISAAC field, we measured the shift with respect to 
the ACS i775-band mosaic for stars and compact sources 
using the imexam task in IRAF. A 2(7-clipped sample 
of reference sources typically consisted of 15-20 objects 
per ISAAC field. The 2(T-clipping was applied to exclude 
sources with strong morphological k-corrections. The dif- 
ference between the shifts implied by individual reference 
sources and the final astrometric solution had a standard 
deviation of less than 0'.'09 in all NIR bands. A map of 
residual shifts for the i^Tg-band mosaic with respect to 
the convolved ACS i775-band mosaic is presented in Fig. 
2. 

First we applied the fractional pixel shift for each 
ISAAC field in each band using the IRAF imshift task 
with a cubic spline interpolation. Next, we summed the 
integer pixel shifted fields applying an identical weight- 
ing scheme as described by FS06 to optimize the S/N for 
point sources, namely: 



I^PIX / NO 

(rmsi.sFWHM)"' 



(1) 



where the weight factor for a given pixel Wpix equals its 
value in the normalized weight map Wnorm, scaled with 
the square of the rms noise measured within apertures of 
1.5FWHM diameter that are placed randomly in empty 
regions of the ISAAC field (see §3.5). Here, uinorm is sim- 
ply the weight map as released by ESO/GOODS normal- 
ized to a maximum of unity. It thus contains information 
on the relative integration time per pixel and accounts 
for bad/hot pixels and pixels with other artifacts that 
were excluded. The weight maps for each field given by 
Eq. 1 were combined into a mosaic that accounts for the 
effective background noise over the entire area. In the 
catalog, we give the weight value for each object mea- 
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surcd on this effective weight map, normahzed to the 
median weight value in the considered band. 

The standard deviation of positional offsets between 
reference sources on the registered WFI images and the 
Ks-hand mosaic were less than O'.'l. 

We chose not to combine the 2 epochs of IRAC observa- 
tions into one mosaic because the 180 degrees difference 
in position angle would lead to a different PSF shape in 
the overlap region than in either of both single epoch ar- 
eas, demanding the use of a different convolving kernel 
over different parts of the field. Instead, we treat each of 
the IRAC epochs independently, providing an empirical 
quality check of the photometry in the overlap region. 
The registration of each of the IRAC images (epoch 1 
and 2) onto the 2x2 blocked Ks-ha,nd image has a po- 
sitional accuracy of better than 0'.'4, as measured from 
offsets between bright star positions on the IRAC and 
iiTg-band images. The positional accuracy for the MIPS 
images is of the order of 0'.'3 rms. We note that minor 
positional offsets between the Kg and IRAC/MIPS im- 
age are solved for by the source fitting algorithm applied 
to IRAC and MIPS photometry (see §4.2.0). 

3.5. Signal to noise and limiting depths 

We analyzed the noise properties of the optical-to- 
24 /nm imaging following the same approach as for 
the FIRES HDFS (LG3) and MS 1054-03 (FS06) data. 
Briefly, the technique uses aperture photometry on 
empty parts of the image to quantify the rms of back- 
ground pixels within the considered aperture size. For 
each convolved ISAAC field in each band, between 200 
and 400 non-overlapping apertures were randomly placed 
at a minimum distance of 5" from the nearest segmen- 
tation pixels in a SExtractor segmentation map. For a 
given aperture size, the distribution of empty aperture 
fluxes is well-fitted by a Gaussian, as illustrated in Fig. 
3(a). We applied a 5a clipping in determining the back- 
ground rms. Panel (b) of Fig. 3 shows that a simple linear 
scaling of the measured background rms a^N) = Na , 
where N = \fA is the linear size of the aperture with area 
A and u is the pixel-to-pixel rms, would lead to underesti- 
mated flux uncertainties. The reason is that correlations 
between neighboring pixels were introduced during the 
reduction and PSF matching. We model the background 
rms as a function of aperture size with a polynomial of 
the form 

^ NHa. + b.N) 

where i refers to the considered band and field, and the 
weight term Wi is derived from the weight map of the re- 
spective field. Fig. 3(b) illustrates the variations in depth 
for the different ISAAC fields, originating from variable 
integration times and observing conditions, and reflected 
in the range of flux uncertainties for objects with simi- 
lar color aperture in the final catalog. For example, the 
upper two curves in the ISAAC Ks panel correspond to 
fields f03 and f04 that had the lowest integration time. 

For the ACS mosaics, we used the same empty aper- 
tures as for the NIR, provided they were within the ACS 
FOV. Every object below the iiT^-band detection thresh- 
old, even though detectable in the ACS imaging, con- 
tributes to the background noise and photometric un- 
certainties of Xg-band detected sources. If we were to 
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Fig. 3. — The background rms derived from the distribution 
of fluxes within empty apertures, (a) Distribution of empty 
aperture fluxes within a 1" , 2" , and 3" aperture diameter on 
the Kg-hanA image of ISAAC field fl5. The distribution is 
well described by a Gaussian with an increasing width for 
increasing aperture size, (b) Background rms as derived from 
flux measurements within empty apertures versus aperture 
size for the ACS bands and the J, H, and Ks ISAAC fields. 
Solid lines represent the functional form from Eq. 2 fit to 
the observed rms noise values. Dashed lines indicate a linear 
extrapolation of the pixel-to-pixel rms. Correlations between 
pixels introduce a stronger than linear scaling with aperture 



restrict our empty aperture analysis to apertures that 
contain neither ii'^-band nor ACS segmentation pixels, 
the background rms estimates for the ACS mosaics would 
decrease by 3 to 9%. In Fig. 3(b), we scaled the back- 



ground rms measured on the ACS and ISAAC images to 
the flux corresponding to AB = 26. 

Uncertainties on the fluxes were assigned based on 
these empirical noise measurements, for an aperture of 
the same area as for the photometric aperture of each 
object. 

Limiting depths of the WFI, IRAC, and MIPS images 
were measured with the same empty aperture method. 
The total limiting AB magnitudes (3cr for point sources, 
i.e., corrected for the flux missed outside the aperture 
used for photometry) are listed in Table 1. 



4. SOURCE DETECTION AND PHOTOMETRY 
4.1. Kg-band Detection 

We aimed to construct a catalog that is especially 
suited to extract stellar mass- limited samples from (e.g., 
van Dokkum et al. 2006). Although the rest-frame NIR, 
probed by IRAC, is a better tracer for stellar mass than 
the rest-frame optical, the downside is its coarser resolu- 
tion, leading to severe confusion. Therefore, we decided 
to detect sources in the observed i^T^-band. 

We used the SExtractor v2.2.2 source extraction soft- 
ware by Bertin & Arnouts (1996) to detect sources with 
at least 1 pixel above a surface brightness threshold of 
^{Ks^ab) — 24.6 mag arcsec"^, corresponding to w 5a of 
the rms background for a typical i^T^-band field. Setting 
the threshold to the same number of ADUs across the 
image instead of adopting a S/N criterion was favored, 
since in the latter case the varying noise properties in the 
i^s-band mosaic would lead to different limiting magni- 
tudes and limiting surface brightnesses from one field to 
the other. We smoothed the detection map with a gaus- 
sian flher of FWHM = 0'.'65, the size of the PSF in the 
detection image. This procedure optimizes the detection 
of point sources. 

The resulting catalog contains 6308 sources, 5687 of 
which have a position on the i^T^-band mosaic where the 
pixel weight is larger than 30% of the median weight. 
We recommend applying such a weight criterion in order 
to construct samples with robust photometry. Running 
SExtractor with identical parameters on the detection 
map multiplied by -1, we obtain a total of 43 spurious 
sources in the area with more than 30% of the median 
weight. Only one of these has S/Nk^ > 5. 

SExtractor flagged 12% of the detected sources as 
blended and/or biased. These sources were treated sep- 
arately in the photometry. 



4.2. Photometry 
Space-based optical and NIR photometry 

We performed the photometry on the convolved ACS 
^435^606*775-2850 and ISAAC JHKs mosaics using SEx- 
tractor in dual image mode, with the i^T^-band mosaic as 
detection map. We derive the color and total aperture 
from the detection image. The same apertures were used 
in each band. We follow L03 and FS06 in defining the 
color aperture based on the Kg-hemd isophotal aperture, 
more precisely on the equivalent circularized isophotal 
diameter diso = 2{Aiso/T^)^''^, where Aiso is the area of 



the isophotal aperture. For isolated sources, we apply 

{APER{ISO), I'.'O < diso < 2'.'0 
APER{1'.'0), diso < I'.'O 
APER{2:.'0), diso > 2'.'0 

(3) 
where APER(ISO) refers to the isophotal aperture 
defined by the surface brightness detection thresh- 
old. Blended sources (indicated with SExtractor flag 
"blended" or "biased" ) were treated separately, 

APER{di,o/s), I'.'O < diso/s < 2'.'0 

APER{1'.'0), diso/s < I'.'O 

APERC^.'O), diso/s > 2'.'0 



APER{COLOR) 



where the reduction factor s for the aperture sizes is in- 
troduced to minimize contamination by close neighbors. 
We adopt the optimal value of s = 1.4 that was deter- 
mined from experimentation by L03 and FS06. 

The motivation for the tailored isophotal apertures de- 
fined in Eq. 3 and Eq. 4 is that it maximizes the S/N 
of the flux measurement. The minimum diameter of I'/O 
corresponds to 1.5 x FWHM of the PSF-matchcd mo- 
saics. The maximum diameter of 2'/0 was adopted to 
avoid flux from neighboring sources and avoid the large 
uncertainties corresponding to large isophotal apertures. 

SExtractor's "MACAUTO" was used to derive the 
total flux of the Kg-haxid detected objects, unless the 
source was blended, in which case the total aperture was 
set to the color aperture: 



APER{TOTAL) 



APER{AUTO), isolated sources 

APER{COLOR), blended sources 

(S) 
Finally, an aperture correction was applied to compute 
the total integrated flux. The correction factor equaled 
the ratio of the flux within an aperture of radius 3" and 
radius rtot — (Aot/77)^/^, based on the curve of growth 
of the PSF constructed from stellar profiles. Here, Atot 
represents the area of the total aperture. The reference 
aperture of 3" radius is sufficiently large since less than 
0.1% of the flux for point sources falls outside this radius. 

Flux uncertainties in both color and total aperture 
were derived from Eq. 2. The quoted uncertainties thus 
take into account both the aperture size used for the flux 
measurement and the limiting depth in the respective 
region of the mosaic. 

We performed a large number of simulations to de- 
termine the accuracy of our measurements of the total 
Kg-hand magnitude. In each simulation, 200 artificial 
sources were added to the X^-band mosaic. They had 
either exponential (Sersic n=l) or de Vaucouleurs (Ser- 
sic n=4) light profiles that were truncated at 5 effec- 
tive radii, and random position angles and ellipticities. 
Source magnitudes and sizes were drawn from the true 
distribution in our catalog in order to span the entire 
range of real galaxy properties. Subsequently, we ran 
SExtractor with an identical configuration as for the real 
catalog, followed by the procedure to compute KI°^q 
as described above. The recovered total magnitudes 
were within 0.1 (0.2) mag from the input value for 89 
(96)% of the sources with exponential light profiles. Since 
our aperture corrections are based on the stellar growth 



curve, our procedure could not account for all flux in the 
extended wings of low surface brightness sources with 
a high concentration (Sersic n=4). The recovered total 
magnitudes were within 0.1 (0.2) mag for 82 (91)% of the 
sources with de Vaucouleurs profiles. The sources where 
the Kl°p^Q underestimate exceeded 0.2 mag generally had 

/ie > 22.5 mag arcsec^^. Taking into account the spread 
in Sersic indices for real sources in our catalog, we cau- 
tion that underestimates of the total i^T^-band magnitude 
by over 0.1, 0.2, and 0.3 mag could occur for roughly 13, 
6, and 4% of the galaxies in our catalog, predominantly 
those with large size and low surface brightness. 

WFI, IRAC and MIPS 24 fim photometry 

In this section, we describe the photometry of i^T^-band 
detected objects in the ground-based optical (WFI) and 
Spitzer IRAC and MIPS 24 ^m imaging of the CDFS. 
For an in depth discussion of the source fitting algorithm 
used, and simulations of its performance, we refer the 
reader to Labbe et al. (in preparation). A short descrip- 
tion with illustration was also presented by Wuyts et al. 
(2007). 

Briefiy, the information on position and extent of 
sources based on the higher resolution Kg-h&nd seg- 
mentation map was used to model the lower resolution 
ground-based optical and 3.6 /im to 24 /j.m images. Each 
source was extracted separately from the Kg-hand image 
and, under the assumption of negligible morphological 
k-corrections, convolved to the WFI, IRAC, respectively 
MIPS resolution. A fit to the WFI/IRAC/MIPS image 
was then made for all sources simultaneously, where the 
fluxes of the objects were left as free parameters. Next, 
we subtracted the modeled light of neighboring objects 
and measured the flux on the cleaned WFI/IRAC/MIPS 
maps within a fixed aperture: 2" for the WFI bands, 3" 
for the IRAC bands and 6" for the MIPS 24 ^m band. 
In order to compute a consistent iCs-WFI/IRAC/MIPS 
color, we measured the source's flux fconv.Ks on a cleaned 
i^s-band image convolved to the WFI/IRAC/MIPS res- 
olution within the same aperture. We then scaled the 
photometry to the same color apertures that were used 
for the space-based optical and NIR photometry, allow- 
ing a straightforward computation of colors over a U^s- 
to-8 ^m wavelength baseline. For the IRAC photometry, 
this means the catalog flux was computed as follows: 



/lRAC,col — /lRAC,3" * 



Ik.. 



col 



conv. Ks ,3" 



(6) 



For the WFI photometry, an identical procedure was fol- 
lowed. An aperture correction based on the growth curve 
of the 24 /im PSF was applied to scale the 24 /im flux 
measurements to the integrated 24 /im flux under the 
assumption that the sources were point-like. 

We note that several other authors developed simi- 
lar algorithms to reduce the effects of confusion (Perez- 
Gonzalez et al. 2005, 2008; G06a). Perez- Gonzalez et 
al. (2005, 2008) compute SExtractor photometry for an 
IRAC-selected catalog that is later merged with shorter 
wavelength photometric catalogs. For the IRAC sources 
with multiple UV/optical/NIR counterparts within 2'.'5, 
the photometry was recomputed using a deconvolution 
method that is similar to that by Labbe et al. (in prepa- 
ration), except for two aspects. The algorithm by Labbe 



et al. (in preparation) takes into account the spatial ex- 
tent of the sources on the reference (Kg-hand) image. 
Moreover, Labbe et al. (in preparation) does not adopt 
the best-fit flux from the source fitting as final photom- 
etry, but rather measures the fiux within an aperture on 
the cleaned image (followed by an aperture correction). 
This allows more robust photometry in cases where the 
object profile varies from the reference to the low reso- 
lution image. A comparison with the multi-wavelength 
photometry by G06a is presented in §7.2. 

Uncertainties in the measured fiuxes in the 
WFI/IRAC/MIPS bands have a contribution from 
the background rms and from the residual contami- 
nation of the subtracted neighbors (Labbe et al., in 
preparation). The former was calculated using the 
empty aperture method described in §3.5. The lat- 
ter was derived by scaling the normalized convolved 
Ks-hand image of each neighboring source with the 
formal la error in its fitted WFI/IRAC/MIPS fiux as 
computed from the covariance matrix produced by the 
least-squares minimization. Subsequently, we performed 
aperture photometry on this image of residual neighbor 
contamination. This gives the confusion error corre- 
sponding to the source flux measured within an aperture 
of the same size. The two contributions were added in 
quadrature to obtain the total error budget. In the 3.6 
/im band, where confusion is most severe, the confusion 
term is typically 5 times smaller than the background 
noise term. For 6% of the sources in our catalog, the 
confusion term dominates the total error budget. 

Here, we follow an empirical approach to validate the 
size of the uncertainties in the IRAC photometry. We 
exploit the overlap region between the 2 independent ob- 
servation epochs of the CDFS with the IRAC instrument. 
The position angle was rotated over 180 degrees, caus- 
ing the PSF to have a different orientation with respect 
to the positions of neighboring sources. In Fig. 4, we 
show the difference between the IRAC magnitude mea- 
sured during epoch 1 and epoch 2. The rms ranges from 
5% in the 4.5 /im band to 10% in the 8.0 /im band for 
sources with an AB magnitude brighter than 22. The 
largest systematic offset was measured for the 3.6 /im 
band, where a zero-point drift of 0.03 mag was measured 
between the 2 epochs. In the inset panels the distribution 

of (/cpochi - /opoch2)/y^err2p^^j^^ + err^^^^^^^ is plotted. 

The distribution is well described by a gaussian. For 
well estimated errors the expected standard deviation of 
the distribution is 1. We adopted a minimum relative 
uncertainty in the flux of 3% to account for zero-point 
variations over the field. This is particularly relevant 
for the 3.6 /im and 4.5 fim band, where the sources are 
detected with a high signal-to-noise. The standard de- 
viation of (/cpochi - .fcpoch2)/y^err2p^^j^^ + err^^^^^^^ in 

these bands is smaller than 1. Adopting a more conser- 
vative minimum relative uncertainty would only decrease 
this value, suggesting that zero-point variations within 
the field are limited to the few percent level. In the less 
sensitive 5.8 /im and 8.0 /im bands, where the minimum 
relative uncertainty is not reached, we find a distribution 

of (/epochi - /cpoch2)/-^err2p^^j^^ _^ err2p^^j^2 ^ith a stan- 
dard deviation of nearly unity, confirming empirically the 
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Fig. 4. — Comparison between IRAC observations from epoch 1 and epoch 2 for A's-band detected sources in the over- 
lap region between the 2 epochs. The large panels show a good correspondence between the 2 independent photometric 
measurements, with a slight zero-point drift of 0.03 mag in the 3.6 /im band. The inset panels shows the distribution of 

(/cpochi — /cpoch2)/\/err^ j^j^ 4- err^ j^j) where a minimum relative error of 3% was assumed to account for relative zero-point 

uncertainties over the field. The standard deviation of the distribution is of order unity, meaning that estimated flux errors 
account well for the empirically determined uncertainties. 



validity of our estimated uncertainties. 

5. REDSHIFTS 

5.1. Spectroscopic redshifts 

The CDFS-GOODS area has been targeted intensively 
by various spectroscopic surveys, listed in Table 3. The 
combined sample of spectroscopic redshifts forms a het- 
erogeneous family of objects, with selection criteria vary- 
ing from pure /-band (VVDS, Le Fevre et al. 2004), 
iiTs-band (Mignoli et al. 2005) or X-ray (Szokoly et al. 
2004) flux limits to various color criteria (e.g., Doherty 
et al. 2005; Wuyts et al. in preparation). It is therefore 
impossible to build a complete Ks-hand selected spec- 
troscopic sample from the data at hand. Rather, we 
aim to provide a list of trustworthy spectroscopic red- 



shifts that are reliably cross-identified with a Ks-hand 
detection in our catalog. To do so, we apply a conser- 
vative quality cut based on the quality flags that come 
with each of the spectroscopic catalogs, and assign the 
redshift to the nearest i^Tg-band selected object within a 
radius of 1'.'2. The quality flags and number of sources 
included in our reliable list of cross-correlated spectro- 
scopic redshifts are summarized in Table 3. We mark 
these sources with a "zsp_qual" flag of 1 in our cata- 
log. For completeness, other spectroscopic redshifts for 
Ks-hand detected objects are also listed in our catalog, 
marked with a "zsp_qual" flag lower than 1 (zsp_qual=0.5 
for F0RS2 quality flag B and VIMOS quality flag B, and 
zsp_qual=0.1 for all others), together with the original 
quality flag from the respective survey. We proceed to 
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Table 3. Spectroscopic hedshifts for Ks-banb detected objects 



Survey 




High quality flags 


Number'^ 


FORS2 (v3.0, Vanzella et al. 


2008) 


A 


357 


K20 (Mignoli et al. 2005) 




1 


265 


VVDS (vl.O, Le Fevre et al. 


2004) 


4,3 


251 


CXO (Szokoly et al. 2004) 




3,2,1 


89 


Norman et al. 2002 




all 


1 


Croom et al. 2001 




all 


21 


van der Wei et al. 2004 




all 


24 


Cristiani et al. 2000 




all 


3 


Strolger et al. 2004 




all 


5 


Daddi et al. 2004 




all 


7 


IMAGES (Ravikumar et al. 


2006) 


1 


100 


LCIRS (Dohcrty ct al. 2005) 




3 


3 


Wuyts et al. (in prep.) 




all 


7 


Kriek et al. (2007) 




all 


2 


Roche et al. (2006) 




all 


4 


Huang et al. (in prep.) 




all 


9 


VIMOS (vl.O, Popesso et al. 


2008) 


A 


329 



''The numbers are non-redundant. For objects targeted during 
multiple surveys, the rcdshift with the highest quality flag was 
adopted. 



use only the 1477 spectroscopic redshifts with zsp_qual 
= 1. 

5.2. Photometric redshifts 

Together with the observed photometry, we release a 
list of photometric redshifts computed with the new pho- 
tometric redshift code EAZY. A full description of the al- 
gorithm and template sets will be presented by Brammer 
et al. (in preparation). Briefly, the program fits a non- 
negative superposition of SED templates to the C/ss-to-S 
/xm photometry, using a template error function that ef- 
fectively downweights the rest-frame UV and rest-frame 
NIR of the templates in the fit. Next, a redshift probabil- 
ity distribution p{z\C, mo) is constructed for each galaxy 
with observed colors C, using the iiT^-band magnitude 
mo as a prior. Finally, the value Zmp of the redshift 
marginalized over the total probability distribution. 



/ Z p(z|C, TTlg) dz 

'"'" Jp{z\C,mo) dz 



(7) 



was adopted as best estimate of the galaxy's redshift. 
la confidence limits were computed by integrating the 
redshift probability distribution from the edges till the 
integrated probability equals 0.317/2. 

The template set constructed by Brammer et al. (in 
preparation) is based on PEGASE models (Fioc & 
Rocca-Volmerange 1997). From a library of ^ 3000 tem- 
plates described by G06a, Brammer et al. (in prepa- 
ration) followed the "non-negative matrix factorization" 
algorithm of Blanton & Roweis (2007) in order to derive 
a set of 5 principal component templates that span the 
colors of galaxies in the semi-analytic model by De Lucia 
& Blaizot (2007). One dusty template with an age of 50 
Myr and Ay = 2.75 was added to account for the exis- 
tence of dustier (and thus redder) galaxies than present 
in the semi-analytic model mock catalog. It is this set of 
6 spectral templates that we fed to EAZY. 

We quantify the accuracy of the photometric red- 
shifts Zphot by a comparison to the spectroscopic red- 
shifts with zsp_qual = 1. Fig. 5(a) shows the corre- 
spondence between Zphot and Zgpcc for all 1067 sources 
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Fig. 5. — Comparison between photometric and spectro- 
scopic redshifts for 1067 T^s-band detected sources with reli- 
able Zapec identification and coverage in all wavelength bands, 
(a) A direct comparison with la confidence intervals derived 
from the redshift probability distribution, (b) The distribu- 
tion of A2/(l + 2:). 2.5% of the sources fall outside the plotted 
range. 



with i^*°AB ^ 24.3 that are covered by all bands and for 
which a reliable spectroscopic redshift is available. The 
error bars indicate the la confidence intervals as derived 
from the redshift probability distribution. We find that 
the spectroscopic redshift lies within the formal la con- 
fidence interval for 66% of the sources, confirming that 
the uncertainties on Zphot have reliably been established. 
Fig. 5(b) presents the distribution of Az/(1 + z), with 



Az 



-2-phot 



which is commonly used to deter- 



mine the accuracy of photometric redshifts. We find a 
median Az/(1 -I- z) of -0.001 and a normalized median 
absolute deviation (equal to the rms for a Gaussian dis- 
tribution) of (Tnmad = 0.032. 3% of the objects with 
spectroscopic redshift have |Az|/(l -I- z) > 5a'NMAD- 
We investigated the quality of Zphot estimates as a func- 
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tion of redshift and total Xg-band magnitude. Variations 
in the median /S.z/(l+z) as function of redshift are small, 
reaching a maximum of 0.05 in the 2.5 < z < 3 interval. 
No obvious dependence on K\°p^^ is observed. The scat- 
ter in Az/(1 + z) is roughly a factor 2.4 larger for the 
z > 1.5 regime (ctnmad — 0.071) than for the z < 1.5 
(o-NMAD = 0.030) regime. 

Considering the 98 spectroscopically confirmed sources 
with a cross-identification within 2'.'0 in the IMs X-ray 
catalog by Giacconi et al. (2002), we find a scatter of 
CNMAD — 0.043. It is reassuring that despite the lack 
of AGN spectrum in our template set, the netto perfor- 
mance of our photometric redshift code for AGN can- 
didates remains good. We do note however that, inde- 
pendent of redshift, the fraction of catastrophic outliers 
(|Az|/(l -I- z) > 5crNMAD) is 3 times larger for the AGN 
candidates than for the total sample of spectroscopically 
confirmed sources. 

6. CATALOG PARAMETERS 

Here we describe the entries of our Ks-ha,nd selected 
FIREWORKS catalog of the GOODS-CDFS. The format 
is similar to the FIRES catalogs of the HDFS (LOS) and 
MS 1054-03 (FS06), making a straightforward combina- 
tion of all three fields possible for the user. The catalog 
can be obtained from the FIREWORKS homepage^. 

• ID- Unique identification number 

• X, y- Pixel position of the object, based on the 
Ks-ha.nd detection map. The pixel scale is 0'.'15 
pixel" ^. 

• RA, DEC- Right ascension and declination coordi- 
nates for equinox J2000.0. 

• [band]_co\i- Flux in microjanskys measured within 
the color aperture (§ 4.2.0). The bandpasses are 
U38,B435,B,V,Vea6,R,'i775,I,Z850,J,H,Ks, [3.6^m], 
[4.5/im], [5.8/im], and [8.0/im]. 

• [band]-Co\ie- Uncertainty in the [band]-Co\i fiux 
measurement, derived from the noise analysis (§ 
3.5). The units are microjanskys. 

• KgJjOti- Total i^T^-band flux in microjanskys, mea- 
sured within the total aperture and scaled by the 
aperture correction (§ 4.2.0). Approximate to- 
tal fluxes in other bandpasses can be computed 
by [band]Jotf = [band]-Colf x {Ks-totf /Kg-colf), 
provided there are no large morphological k- 
corrections. 

• iiTs.totfe- Uncertainty associated with i^s_totf, also 
in microjanskys. 

• [24/xm]_totf- Total MIPS 24 ^m-band fiux in mi- 
crojanskys, measured within a 6" diameter circular 
aperture and then aperture corrected (§ 4.2.0). 

• [band\w- Effective weight in the bandpass [band], 
normalized to the median effective weight of all 
sources in that band. We recommend applying a 
conservative weight criterion [band]w > 0.3 to con- 
struct samples with robust photometry. 

^ http://www.strw.lcidcnuniv.nl/fircworks 



• ap_col- Aperture diameter in arcsec within which 
[6aK(i]_colf was measured. In cases where the color 
aperture was the isophotal aperture defined by the 
surface brightness threshold of h{Ks,ab) = 24.6 
mag arcsec"^, ap_col is the diameter in arcsec of a 
circular aperture with equal area. 

• ap_tot- Aperture diameter in arcsec used for mea- 
suring the total Xg-band flux. When the isophotal 
or SExtractor's "MAG_AUTO" aperture was used, 
this entry contains the equivalent circularized di- 
ameter corresponding to that aperture. 

• f_deblendl- Flag equal to 1 when the source was 
deblended somewhere in the process (SExtractor's 
"blend"). 

• f_deblcnd2- Flag equal to 1 when the photometry 
is affected by a neighboring source (SExtractor's 
"bias"). 

• Kr50- Half light radius in arcsec, measured on the 
Ks-ha.nd image (SExtractor's flux_radius scaled to 
arcsec) . 

• Kcps- EUipticity of the isophotal area, measured 
on the Ks-hand image. 

• Kposang- Position angle of the isophotal area, 
measured on the Ks-ha,nd image. 

• zph_best- Best estimate of the photometric redshift 

(§ 5.2). 

• zphJow, zph_high- Lower and upper edge of the 
68% confidence interval around zph_best. 

• zsp- Spectroscopic redshift (set to -99 when no 
spectroscopic information is available). 

• zsp_qual- Quality flag from to 1 assigned to the 
spectroscopic redshift. Only zsp_qual=l entries are 
considered reliable. 

• zsp_source- Spectroscopic survey from which zsp 
was taken (Tab. 3). 

• zsp_quaLorig- Original quality flag for zsp from the 
respective spectroscopic survey. 

• XID- Identification number from the IMs X-ray 
catalog by Giacconi et al. (2002), set to -99 when 
no cross-identification within 2" was found. Note 
that we accounted for the ~ 1'.'3 systematic offset 
in the Giacconi et al. (2002) X-ray centroids, as 
pointed out before by Roche et al. (2003). 

7. COMPARISON TO THE GOODS-MUSIC CATALOG 

7.1. Differences in data and strategy 

Recently, G06a presented a multicolor catalog for the 
GOODS-CDFS field, referred to as the GOODS-MUSIC 
catalog. The clustering evolution of distant red galax- 
ies was quantified based on this sample (Grazian et al. 
2006b), as was the contribution of various color-selected 
samples of distant galaxies to the stellar mass density 
(Grazian et al. 2007). Despite the overlap in public data 
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used to compile the GOODS-MUSIC and our catalog, 
there are a number of marked differences. 

First, our catalog is purely Ks-ha,nd selected. Since 
the Ks-hand magnitude is a better proxy for stellar mass 
than the optical magnitude, this makes it ideally suited 
to extract mass-limited samples from. The GOODS- 
MUSIC sample on the other hand is to first order zs5o- 
band selected (at the ACS resolution), with an addition 
of the remaining ii'^-band sources that are detected in a 
map with masked 2;85o-band detections. Although valu- 
able in its own respect, this makes it less trivial to un- 
derstand the completeness of the sample. 

Second, we based our catalog on the ESO/GOODS 
data release vl.5, consisting of 3 extra ISAAC pointings 
in J and Kg, and 7 more in H with respect to the vl.O 
release used by G06a. 

Third, apart from the U^s photometry, we also added 
photometry based on images in the other WFI pass- 
bands: BVRI. Despite the shallower depth than the 
ACS imaging (see Table 1), this addition is useful for, 
e.g., estimating photometric redshifts, since the WFI 
passbands are nicely centered in between the ACS pass- 
bands. Consequently, the spectral energy distributions 
of sources in the catalog are better sampled. 

Finally, we include MIPS 24 /xm measurements, en- 
abling us to constrain the total IR luminosity of the Kg- 
band selected galaxies. Before doing so, we compare the 
photometry in common between both catalogs, and the 
photometric redshifts derived from it. 

7.2. Comparing photometry 

We cross-correlated the two catalogs using a search ra- 
dius of 1'.'2 and in Fig. 6 present a comparison of the [/as- 
to-8.0 /xm photometry for the bands in common between 
them. Differences between the total magnitudes are plot- 
ted for objects with S/N > 10 in the i^g-band and the 
band under consideration. Objects that are marked by 
SExtractor as blended in the Kg-hand, are indicated with 
empty symbols. 

Our Uss magnitudes for non-blended sources are 
fainter than the GOODS-MUSIC magnitudes by 0.2 mag 
in the median, with a scatter of 0.4 mag. This offset can 
be attributed to our use of the re-reduced WFI data for 
which the GaBoDS consortium also supplied an improved 
zeropoint determination. 

The overall correspondence in the B-to-Kg bands is 
good, and offsets can be well understood from the dif- 
ferences in the applied photometric method. We mea- 
sure a typical median offset for non-blended sources 
in the optical and NIR bands of rnagtot, fireworks ~ 
"lagtot, MUSIC == -0.06, and a scatter of ctnmad < 
0.2. G06a based their total magnitudes on SExtractor's 
"MAG_AUTO" parameter for the zgso-band detections 
and on the "MAG.BEST" for the remaining JsT^-band de- 
tections that were not detected in the zsso-band. G06a 
did not apply an aperture correction based on the stellar 
growth curve to correct for the flux lost because it fell 
outside the "MAG_AUTO" or "MAG.BEST" aperture. 
The lack of aperture correction explains at least part of 
the systematic offset. Sources marked as blended in our 
Kg-hand detection map typically are brighter by 0.2 - 
0.4 mag in the MUSIC catalog. This can be explained 
by the contamination from neighboring sources within 
the "MAG_AUTO" aperture, which we avoid by using 



the isophotal aperture in combination with an aperture 
correction for blended sources. However, our measure of 
the total magnitude in the presence of crowding is also 
not perfect in that flux from the source that is projected 
on top of a blended source may be lost by our technique. 
For the IRAC photometry, the discrepancies are larger, 
ranging from 0.16 mag in the 4.5 /im band to 0.42 mag 
at 8.0 fiiai. The observed offsets in IRAC photometry 
can largely be attributed to the use of an early version 
of the IRAC PSF by G06a and a bug in the normal- 
ization of the smoothing kernel for the IRAC data by 
G06a (Fontana & Grazian, private communication). Af- 
ter removing these two effects, our measurements of the 
total IRAC magnitudes are still brighter by 0.06 to 0.1 
mag. This offset is similar to that for the optical and 
NIR bands, and can be attributed to light missed by the 
"MAG.BEST" magnitude, as G06a aheady caution. We 
stress that the offset in total IRAC magnitudes not only 
affects its derived properties such as stellar mass, but also 
the optical-to-MIR and NIR-to-MIR colors. For exam- 
ple, our 2:850 -[3. 6^m], zgso- [4.5/xm], zsso- [5.8^m], and 
Z850- [S.O^m] colors are redder than the GOODS-MUSIC 
colors by 0.23, 0.11, 0.17, and 0.37 mag in the median 
respectively. Similarly, our Kg — [3.6/im], Kg — [4.5/im], 
Kg — [5.8/im], and Kg — [8.0/im] colors are redder in the 
median by 0.29, 0.16, 0.22, and 0.42 respectively. The 
scatter in the color differences with respect to GOODS- 
MUSIC typically amounts to 1.5 times the size of the 
median offset. 

7.3. Comparing photometric redshifts 

Finally, we compare the photometric redshifts pre- 
sented in §5.2 with those derived by G06a. The numbers 
quoted in §5.2 and by G06a cannot directly be compared 
since new spectroscopic redshifts were added, and objects 
that showed evidence for the presence of an AGN in their 
optical spectrum were rejected from the GOODS-MUSIC 
photometric redshift analysis. Nevertheless, when com- 
paring the performance of the Zphot estimates for a set 
of 659 non-AGN with reliable Zgpec and coverage in all 
bands in both catalogs, we find a scatter in Az/{l + z) for 
GOODS-MUSIC (fJNMAD = 0.039) that is similar to that 
for our best estimates (ctnmad = 0.033). The median 
Az/(1 + z) is -0.009 and -0.005 for the GOODS-MUSIC 
and FIREWORKS Zphot estimates respectively. 

Comparing the Zphot estimates for all objects that 
are well exposed in all bands in both catalogs, ir- 
respective of a spectroscopic confirmation, we find a 

median ^Pt.°t music -^phot.FIREWORKS Qf Q ^^^ ^ g^^^^^gj. 

-L + S^phot, FIREWORKS 

CTNMAD — 0.052. The largest systematic offsets occur 
in the 1.5 < z < 2 and 2.5 < z < 3 intervals where 
median( ^p^ot music -^phot.FiREwoRKs ) _ q Q4 ^^^ _q Q4 ^^_ 

^ i4-2phot, FIREWORKS ' 

spectively. We find no systematic dependence on K\°\^. 
We conclude that the photometric differences be- 
tween both catalogs can be understood from the applied 
method. Our Zphot estimates are in excellent agreement 
with the available spectroscopic samples and generally 
show good agreement with the Zphot estimates by G06a. 
We now proceed to exploit our catalog to analyze the 
colors and total IR energy output of distant galaxies. 

8. TOTAL IR PROPERTIES OF DISTANT Xs-SELECTED 
GALAXIES 
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Fig. 6. — A direct comparison of total magnitudes for sources with S/N > 10 in the C/as-to-S.O pm bandpasses in common 
between GOODS-MUSIC and our FIREWORKS catalog. Sources that are blended in the iC's-band image are plotted as empty 
symbols. On the right side of each panel, a histogram shows the distribution of offsets. Our U^s photometry, based on the re- 
reduced WFI image by the GaBoDS consortium, is systematically fainter by 0.2 mag. We find an overall good correspondence in 
the ACS optical and ISAAC NIR bands, with offsets of roughly 6% due to aperture corrections. In the IRAC bands, photometric 
offsets amount up to 0.4 mag. 



With the catalog at hand, we aim to answer the fol- 
lowing simple questions: Which iiTs-selected (S'/A'^i^:^ > 5, 
-^s°AB ^ 24.3) galaxies at 1.5 < z < 2.5 have the bright- 
est total IR luminosities (Lir = L{8 — 1000 /^m)), and 
which contribute most to the integrated total IR lu- 
minosity? Specifically we will address whether the to- 
tal IR luminosity is dominated by red or blue galax- 
ies, with the color defined in the rest-frame UV, opti- 
cal or NIR wavelength regime. We restrict our analy- 
sis to the 1.5 < z < 2.5 interval, since at those red- 
shifts the observed 24 /xm probes the rest-frame mid-IR 
(7 /im< Arest ^ 10 /xm), which broadly correlates with 
the total IR luminosity (e.g., Spinoglio et al. 1995; Chary 



& Elbaz 2001; Dale & Helou 2002). 

8.1. Observed 24 IJ-m flux as function of observed colors 

We approach the questions raised above by first study- 
ing the correlation between purely observational proper- 
ties: the 24 fim fiux as proxy for IR luminosity and the 
observed i?435 — Veoe, J ~ Ks, and Ks — [4.5/im] colors as 
proxy for the rest-frame UV, optical and optical-to-NIR 
color respectively. Unless the redshift dependence of the 
conversion from 24 /im to total IR luminosity and of the 
conversion from observed to rest-frame colors are con- 
spiring, any trend in the directly observable properties 
should be a signpost for correlations in the rest-frame 
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properties, whose derivation involves significant system- 
atic uncertainties. 

The total number of i^Tg-selected galaxies at 1.5 < z < 
2.5 with Kl°lj^ < 24.3 is 961. 45% of these are indi- 
vidually detected at the 3cr level ([24/im]tot > 11 /^Jy). 
Since we aim to investigate what color galaxies dominate 
the IR emission, we divide our galaxies in bins of simi- 
lar color. Each bin contains an equal number of sources 
(160). Consequently, the bin widths are not equal. To 
start, we leave the origin of the 24 /xm emission (dust 
heated by AGN or star formation) as an open question. 
We note however that when excluding X-ray detected 
sources each bin would contain 154 objects, and apply- 
ing such a selection would not affect the results of our 
analysis. We perform the stacking per color bin in two 
ways: by averaging the 24 /xm flux measurements of all 
sources in the bin (some of which will be negative due to 
noise), and by taking the median of the 24 /im flux mea- 
surements of all sources in the bin. The mean has a con- 
tribution from all the galaxies in the color bin, and hence 
directly measures the contribution to the integrated 24 
/xm emission. The median has the advantage that it is 
more robust against a small number of bright sources in 
the color bin, and is therefore lower. 

In Fig. 7, the stacked 24 /im flux densities are plotted 
versus the observed S435 — Veoe; J— Kg, and i^T^ — [4.5/im] 
color. The error bars on the mean flux measurement indi- 
cate the errors in the mean (cr([24]tot)/v^), whereas the 
error bars on the median flux measurement are computed 
as crNMAD([24]tot)/V^- Furthermore, the hght grey and 
dark grey polygon show the range containing respectively 
the central 68% and 50% of the binned galaxies. Each 
color bin contains galaxies with a large spread in 24 /im 
fluxes. In most bins, at least 20% of the galaxies have an 
individual signal-to-noise ratio S/N < 1 at 24 /xm. 

Fig. 7(a) shows that the galaxies in the bluest i?435 — 
14)06 bins are the faintest 24 /xm sources. However, the 
stacked [24/xm]tot flux is not uniformly increasing over 
the whole observed optical color range. Considering col- 
ors measured at longer wavelengths, we find a highly 
significant increase in the stacked [24/im]tot flux over 
the entire J - Kg (Fig. 7(b)) and Kg - [4.5/xm]tot (Fig. 
7(c)) color range. The trend is strongest in the ob- 
served Kg — [4.5/xm]tot color, where we find an increase 
in [24/xm]tot of a factor 26 in the mean and 2 orders of 
magnitude in the median over a color range of ~ 1.3 mag. 
Since the bins contain an equal number of objects, it is 
trivial to see that not only the reddest galaxies in J — Kg 
and Kg — [4.5/xm] are brightest at 24 /xm, they also con- 
tribute the most to the total 24 /xm emission integrated 
over all -ftTs-selected galaxies at 1.5 < z < 2.5. 

8.2. Total IR luminosity as function of rest-frame colors 

Although the trend of more 24 /xm emission for galax- 
ies with a redder observed color is highly significant for 
J — Kg and Kg — [4.5/xm], it could still be contaminated 
or, alternatively, driven by redshift dependencies within 
the 1.5 < z < 2.5 redshift interval under consideration. 
Now, we will attempt to remove possible redshift depen- 
dencies by converting both axes to a rest-frame equiv- 
alent. Moreover, instead of converting the measured 
flux density at 24 /im to a rest-frame fiux density at 
24 /xm/(l -\- z), we use it as a probe to determine the 
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Fig. 7. — Stacked 24 /xm flux densities as function of observed 
B435 — Veoe, J — Ks, and Ks — [4.5/ini] color for galaxies at 
1.5 < 2 < 2.5 with S/N > 5 in the T^s-band (corresponding 
to A'a°AB ^ 24.3). Filled circles represent the median for 
each equal-number bin. Empty boxes represent the mean 
flux. Light-grey and dark-grey polygons indicate respectively 
the central 68% and 50% of the distribution within each bin. 
We find a trend of increasing [24/im]tot with redder observed- 
frame color that gets progressively stronger as we consider 
colors measured at longer wavelengths. 
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total IR luminosity Ljr = i(8 — 1000 /im). Since this 
conversion assumes that the 24 /xm emission originates 
from dust heated by star formation, we further reject 
all X-ray detections from our sample to rule out rel- 
atively unobscured AGN candidates. We caution that 
among the brightest 24 /zm sources in our sample, there 
may still be cases where hot dust emission from an ob- 
scured AGN contributes significantly to the 24 /im emis- 
sion and causes an overestimate of Lir. Detailed studies 
of such mid-IR excess sources are presented by Daddi et 
al. (2007) and Papovich et al. (2007). Papovich et al. 
(2007) find that for sources with [24^m]tot > 250 fiJy 
the LiR estimates derived from the 24 /xm flux density 
alone are too large by factors 2-10. We stress however 
that we study a deep Kg- and not 24 /im-selected sam- 
ple. The fraction of sources with [24/xm]tot > 250 /iJy 
per color bin never exceeds 7%. Therefore, at least the 
median stacked properties are relatively robust against a 
possible contamination by obscured AGN. 

In the following, we first describe the derivation of rest- 
frame UV to NIR colors. Next, we explain the method to 
estimate the total IR luminosity. Finally, we repeat the 
stacking analysis using the derived rest-frame properties. 

Stellar mass, UV slope, and rest-frame colors 

For each of the galaxies in our sample, wc modeled the 
spectral energy distribution (SED) using the stellar pop- 
ulation synthesis code by Bruzual & Chariot (2003) . We 
used an identical approach as Wuyts et al. (2007), as- 
suming a Salpeter IMF and solar metallicity, and fitting 
three star formation histories: a single stellar population 
without dust, an exponentially declining model with e- 
folding time of 300 Myr and allowed dust attenuation 
in the range Ay = — 4, and a constant star forma- 
tion model with the same freedom in attenuation. We 
subsequently scaled the stellar masses derived from the 
best-fitting template to a Kroupa IMF by dividing the 
stellar masses for a Salpeter IMF by a factor 1.6. 

We characterize the rest-frame UV part of each SED 
by fitting the functional form Fx ^ X^ to the best- 
fitting template, using the rest-frame UV bins defined 
by Calzetti, Kinney,& Storchi-Bcrgmann (1994). The 
robustness of this technique is discussed by van Dokkum 
et al. (2006). 

The rest-frame {U — F)rcst and {V — J)rcst colors were 
determined by interpolation between the directly ob- 
served bands using templates as a guide. For an in depth 
discussion of the algorithm, we refer the reader to Rud- 
nick et al. (2001; 2003). We used an IDL implementation 
of the algorithm by Taylor et al. (in preparation) dubbed 
"InterRest" . 

Converting 24 l^-m flux to total IR luminosity 

At redshifts 1.5 < z < 2.5, the 24 ^m fiuxes trace 
the rest-frame 7.7 /im emission from polycyclic aromatic 
hydrocarbons (PAHs). To convert this MIR emission to 
a total IR luminosity Lir, = L{8 — 1000 /im), we use 
the infrared spectral energy distributions of star-forming 
galaxies provided by Dale & Helou (2002). The tem- 
plate set allows us to quantify the IR/MIR flux ratio for 
different heating levels of the interstellar environment, 
parameterized by dM{U) ~ U~°'dU where M(U) rep- 
resents the dust mass heated by an intensity U of the 
interstellar radiation flcld. 



We computed the total infrared luminosity £ir,q for 
each object for all Dale & Helou (2002) templates within 
the reasonable range from a = \ for active galaxies to 
a = 2.5 for quiescent galaxies. The mean of the result- 
ing log(LiR,Q=i^i.o625 2.5) was adopted as best estimate 

for the IR luminosity, and the ±0.45 dex variation from 
-^IR.Q=2.5 to iiR,Q=i was taken as a measure for the sys- 
tematic uncertainty in the conversion. 

Apart from the random photometric error and system- 
atic template uncertainty, uncertainties in the photomet- 
ric redshift contribute to the total error budget. For each 
galaxy, we calculated the spread in Lir caused by vari- 
ations of Zphot within the 68% confldence interval. Al- 
though the uncertainty in photometric redshift is partly 
random (propagating from photometric uncertainties in 
the SED), we treat it as purely systematic, originating 
from template mismatches. This means the error bars 
related to 2;phot on the stacked Lir measurements do not 
scale with l/vN. Instead, they range from the stacked 
Lir based on the lowest LiR^i„dividuai estimates allowed 
for each object within its 2:phot uncertainty, to the stacked 
Lir based on the maximum Lir_ individual allowed for each 
object. 

Lir versus rest-Jrame color 

Having determined the rest-frame colors and total IR 
luminosities of iiTs-selccted galaxies without X-ray de- 
tection at 1.5 < z < 2.5, we now proceed to investi- 
gate which Xg-selected galaxies contribute most to the 
IR emission. Again, we investigate the average proper- 
ties of galaxies binned by color to enhance the robustness 
of our results. In Fig. 8, we plot the mean and median 
total IR luminosities of our in color bins divided sample 
versus the rest-frame UV slope f3, the rest-frame opti- 
cal color {U — V)rost, and the rest-frame optical-to-NIR 
color {V — J)rost- The black error bars indicate the er- 
ror in the mean, respectively median. With dotted error 
bars, we show the systematic variations allowed within 
the photometric redshift uncertainties. Finally, the error 
bar in the bottom right corner represents the range from 
quiescent to active galaxy templates by Dale & Helou 
(2002). Clearly, systematic uncertainties are dominating 
the error budget in this analysis. 

As in Fig. 7, we flnd a large range in IR properties 
within each bin, illustrated by the light grey and dark 
grey polygons that mark respectively the central 68% 
and 50% of the distribution of Lir of individual objects. 
Nevertheless, a general trend is visible of redder colors 
corresponding to larger IR luminosities. No matter which 
part of the spectral energy distribution is used to de- 
flne red or blue galaxies, the redder half of the galaxies 
in our sample have stacked IR properties in the LIRG 
(Lm. — 10^^ — 10"^^ Lq) regime. The trend seems to flat- 
ten at the reddest UV slopes, and the contribution to the 
total IR luminosity even drops for the reddest {U — V)rcst 
bin. In {V — J)rost on the other hand, the increase in 
Lir with reddening color continues over the entire color 
range, reaching ULIRG luminosities (Lir > 10^^ Lq) for 
23% of the galaxies in the reddest bin. 

Considering the IR luminosities of individual objects 
over the entire color range, we flnd that ULIRGs make 
up 7% of our Xg-selected sample with Kl°l^ < 24.3. The 
fraction of ULIRGs increases to 35% when only consider- 
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Fig. 8. — Stacked total IR luminosities as function of rest-frame 
UV slope, (U - y)rest, and {V - J)rest for non-AGN at 1.5 < 
z < 2.5 with S/N > 5 in the K^-band (ii'*°AB £ 24.3). Filled 
circles represent the median for each equal-number bin. Empty 
boxes represent the mean flux. Light-grey and dark-grey polygons 
indicate respectively the central 68% and 50% of the distribution 
within each bin. The systematic uncertainty induced by template 
uncertainties in the conversion to Lij^ is indicated in the bottom 
right corner. The dotted error bars indicate the variation in the 
mean and median Ljj^ by systematic variations in Zpijot- We find 
a trend of increasing Liji with redder rest-frame color. Since the 
bins contain equal numbers of objects, this also means that red 
galaxies contribute most to the integrated IR emission of distant 
if .-selected galaxies. 



ing massive (M > 10^^ Af©) galaxies. A similar fraction 
was found by Daddi et al. (2007). 

The IR/MIR conversion factor varies by nearly an or- 
der of magnitude between the use of quiescent (a = 2.5) 
or active (a = 1) galaxy templates. For each of the 
galaxies, we used the mean of the logarithm of all con- 
version factors derived from the 1 < a < 2.5 templates. 
Other authors (e.g. Papovich et al. 2006; Papovich et 
al. 2007) adopted a conversion to total IR luminosity 
in which the rest-frame IR luminosity density translates 
uniquely to a single SED template (i.e., a lower a was 
used for brighter galaxies). Applying the conversion by 
Papovich et al. (2007) to each of the galaxies in our 
sample and repeating the stacking procedure, gives the 
following results. The mean Ljr for each color bin in- 
creases by typically ~ 0.1 dex with respect to the value 
obtained with our universal (flux-independent) conver- 
sion. The median Lir for each color bin instead decreases 
by typically ^ 0.15 dex. All the described correlations 
between Ljr and rest-frame color remain intact. We thus 
confirm that our results are robust to alternative conver- 
sions from 24 ^m flux to total IR luminosity that are 
commonly used in the literature. 

We conclude that amongst distant Xg-selected galax- 
ies that show no sign of relatively unobscured AGN at 
X-ray wavelengths, the redder galaxies have on average 
larger total IR luminosities. Given that each bin in Fig. 8 
contains an equal number of objects, it is also clear that 
red galaxies in our sample contribute more to the inte- 
grated total IR luminosity than blue galaxies. We argue 
that this trend cannot be explained by systematic errors, 
and tested that this conclusion is robust against the pre- 
cise choice of rcdshift interval by varying the lower edge 
between redshift 1 and 2 and the upper edge between 2 
and 3. Likewise, we verified that none of our conclusions 
critically depend on the number of color bins. Finally, 
we note that, although X-ray detections were excluded 
from our analysis to validate the use of IR templates for 
starforming galaxies, the results remain nearly unaffected 
when we treat them as normal galaxies. 

Dividing the sample in two equal-number bins accord- 
ing to /3, {U — I^)rost7 and {V — J)rcst, we find that the in- 
tegrated total IR luminosity of the red half is larger than 
that of the blue half by a factor 6.7 ± 1.6, 3.6 ± 0.8, and 
5.6 ± 1.2 respectively. Imposing a brighter cut in the Ks- 
band magnitude (Fig. 9) weakens the correlation of Lir 
with UV slope and (F-J)rost- Adopting a ^^"Xb < 22.86 
cut, as is the case for the NIR-selected sample studied by 
Reddy et al. (2006), causes the stacked Lir of the bluest 
(U — T^)rost colors to increase while the stacked Ljr of 
the reddest {U — V)rcst colors drops. The overall fraction 
of galaxies that is not detected at the 3a level in the 24 
/im band is lower by a factor 2 in the Kl°^^ < 22.86 
sample than in our Kl°^^ < 24.3 sample. However, at 
the reddest optical colors {{U — V)rcst > 1-3) the fraction 
of galaxies that is not detected at the 3a level in the 24 
/im band remains the same as we impose the brighter 
Xg-band limit. We find that, for this i^Tg-bright sample, 
the ratio of the integrated Ljr of the red and the blue 
half of the galaxies in f3, {U — T^)rost, and {V — J)rest 
amounts to a factor of 2.1 ± 0.5, 1.2 ± 1.1, and 2.1 ± 0.5 
respectively. 

It is tempting to elaborate on the physical interprcta- 
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Fig. 9. — Stacked total IR luminosities as function of rest- 
frame UV slope, {U — y)rost, and (V — J)rost for non-AGN at 
1.5 < 2 < 2.5 with Kl°lj^ < 22.86 (ii'*°vcga £ 21). Filled circles 
represent the median for each equal-number bin. Empty boxes rep- 
resent the mean flux. Light-grey and dark-grey polygons indicate 
respectively the central 68% and 50% of the distribution within 
each bin. The systematic uncertainty induced by template uncer- 
tainties in the conversion to Ljpi is indicated in the bottom right 
corner. The dotted error bars indicate the variation in the mean 
and median Ljpj^ by systematic variations in z^phot- The increase in 
Li^ with redder rest-frame color is less pronounced than for the 



tion in terms of star formation rate (SFR), age, and dust 
content of the galaxies in our sample implied by the pre- 
sented results. Colors in different wavelength regimes are 
to a greater or lesser extent determined by these physical 
parameters. The UV slope and (T^— J)rcst color are both 
sensitive tracers of dust attenuation (Meurer et al. 1999 
and Wuyts et al. 2007 respectively). In comparison to 
the (V^— J)rest color, the {U—V)rest color is more sensitive 
to stellar age and to lesser extent affected by dust (Wuyts 
et al. 2007). The fact that in the reddest (U-V)^cst bin 
the total IR luminosity drops again might therefore in- 
dicate an increasing contribution from low Lir galaxies 
with little dust-obscured star formation. In combination 
with the fact that rest-frame optically selected galaxies 
often have faint UV luminosities and thus little unob- 
scured star formation (Forster Schreiber et al. 2004), this 
suggests that part of the galaxies making up the reddest 
{U—V)rest bin have a low overall SFR (unobscured -I- ob- 
scured) . This is consistent with the spectral evidence for 
galaxies with quenched star formation at 2; ~ 2 (Kriek et 
al. 2006). A similar conclusion was drawn by Reddy et 
al. (2006), who found for a sample of galaxies at similar 
redshifts selected by optical and NIR color criteria that 
the IR luminosity of 24 iim detected sources increased 
toward redder observed z — K, but that at the reddest 
z ~ K colors a population without 24 /xm detection ex- 
ists that satisfies the Distant Red Galaxy (Franx et al. 
2003) and/or BzK/PE (Daddi et al. 2004) color crite- 
ria. Based on X-ray stacking in the GOODS-North field, 
although also probing only to i^*°AB < 22.86, Reddy et 
al. (2005) found a similar turnover in inferred SFR at 
z ~ K > 3. Although the observed z — K color is red- 
shift dependent and at z ~ 2 spans a somewhat broader 
wavelength range than [U — F)rcst, both colors probe 
the age-sensitive Balmer/4000A break and the observed 
turnovers therefore likely share the same origin. 

An in depth analysis of the mix between dust-obscured 
starforming systems and evolved red galaxies would re- 
quire a careful SED modeling, estimating the SFR based 
on different wavelength tracers from X-ray over UV to 
the infrared, and a treatment of each object on an indi- 
vidual basis to assess their relative contribution. Such 
a study is clearly beyond the scope of this paper, and 
will be presented by Labbe et al. (in preparation) based 
on the combined sample of galaxies in the CDFS (this 
paper), MS 1054-03 (FS06) and the HDFS (L03). 

Redshift dependence 

We measure a median rest-frame color evolution for 
galaxies in our KI°^q < 24.3 sample from z = 2.5 to 
z = 1.5 of 0.48, -0.05, and 0.48 for /?, {U-V),cst, and {V- 
J)rest respectively We note that going to lower redshifts, 
the evolution in (3 and (T^— J)rost reverses. The observed 
evolution amounts to 26, 5, and 40% of the scatter in the 
respective rest-frame color at any given redshift in the 
1.5 < z < 2.5 range. Consequently, not every color bin in 
Fig. 8 and Fig. 9 is composed of galaxies with exactly the 
same redshift distribution. This is particularly the case 
for /3 and (T^ — J)rcst, where the measured color evolution 
with redshift is largest. A dependence on redshift can 
potentially affect the observed trend of increasing iiR 
toward redder colors, if ijR evolves within the 1.5 < 
z < 2.5 interval. Binning galaxies by redshift, we find 
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a decrease in the stacked Ljr from z — 2.5 to z — 1.5 
by ^ 0.4 dex. We conclude that redshift dependencies 
within the 1.5 < z < 2.5 interval are not driving, and 
if anything weaken, the observed increase of Lir toward 
redder colors. 

Mass dependence 

We now investigate the role of stellar mass in the ob- 
served trend of red galaxies dominating the total IR 
emission, van Dokkum ct al. (2006) showed that the 
high-mass end of the high redshift galaxy population is 
dominated by galaxies with red {U — X^)rost and large 
f3. Considering the optical-to-NIR colors of galaxies in 
our sample, a similar lack of massive galaxies with blue 
{V — J) rest is observed. 

In order to investigate whether the Lir, - color rela- 
tion discussed in §8.2.0 merely reflects a correlation be- 
tween stellar mass and Lir, we now divide out the mass 
dependence. Fig. 10 shows the stacked Lir/M, mea- 
surements as a function of color for the -ft'*°AB < 24.3 
sample and the -ft'*°AB < 22.86 subsample. The y-axes 
span 4 orders of magnitude, as in Fig. 8 and Fig. 9, al- 
lowing an easy visual comparison of the strengths of the 
trends. A weaker trend, by up to an order of magnitude, 
is observed between Lir/M* and rest-frame color than 
between Ljr and rest-frame color (Fig. 8 and Fig. 9). 
For the Xg-bright subsample {Kl°l^ < 22.86), the ob- 
served LiR versus /3 and Lir versus {V — J)rost relation 
can even be entirely attributed to a dependence of color 
on stellar mass. Finally, the drop in stacked Ljr/M* at 
the reddest {U — T^)rost colors is stronger than the drop 
in stacked Lir. The population of optically red galaxies 
with only little IR emission, suggested by the analysis in 
§8.2.0, therefore concerns galaxies at the high-mass end 
of the mass function. 

9. SUMMARY 

We present a Xg-band selected catalog for the 
GOODS-CDFS, dubbed FIREWORKS, containing con- 
sistent photometry in the U^s, S435, B, V, Veoe, R, 1775, 
/, Z850, J, H, Kg, [3.6^m], [A.5fim], [5.8^m], [S.Ofim], 
and [24/im] bands. Together with the photometry, we 
release a list of photometric redshifts with a scatter in 
A2;/(l + z) of 0.032, a cross-correlation with all available 
spectroscopic redshifts to date, and a cross-correlation 
with the IMs X-ray catalog by Giacconi et al. (2002). 
After a description of the catalog construction, we dis- 
cuss the differences with the GOODS-MUSIC 'zsso + 
iiTs '-selected catalog by G06a. 

A previous version of the catalog, lacking the WFI 
UssBVRI photometry, has been used to estimate stel- 
lar mass densities (Rudnick et al. 2006), construct lumi- 
nosity functions (Marchesini et al. 2006) and study the 
predominance of red galaxies at the high mass end (van 
Dokkum et al. 2006). In this paper, we exploit the full 
catalog to answer the following question: Which distant 
iiTs-band selected galaxies are brightest and contribute 



most to the total IR luminosity? 

First, we compared the stacked 24 fim fluxes of galaxies 
atl.5<z<2.5 with -ftr' ab < 24.3 split in observed color 
bins. Overall, a large spread in IR properties is found in 
each color bin. Nevertheless, stacking the fluxes within 
each bin reveals a clear trend with color. Both in the 
observed B435 — Vgoe, J — Kg, and Kg — [4.5^m] colors, 
the lowest mean and median [24/im] tot fluxes are found 
for the bluest color bin. In J — Kg and Kg — [4.5/im], the 
emission at 24 /im continues to rise toward redder colors. 

Second, we use our photometric redshifts to convert the 
observed spectral energy distributions to rest-frame col- 
ors and translate the observed 24 /im flux to the total IR 
luminosity Ljr = L(8 — 1000 /xm). In this procedure, all 
relatively unobscured AGN candidates, selected by their 
X-ray detection, were rejected from the sample. Remov- 
ing the redshift dependence and extrapolating from MIR 
to total IR comes at the cost of systematic uncertain- 
ties. We carefully measured the systematic contribution 
to the total error budget from uncertainties in Zphot and 
from our lack of knowledge about which IR template 
SED matches best the spectral shape of the objects in 
our sample. Doing so, we flnd a continuous increase in 
Lir with {V — J)rost- An increasing Lir is also mea- 
sured with UV slope /3, flattening at the largest j3. The 
rising trend of the stacked Lir luminosity toward red- 
der {U — y)rest seems to reverse in the reddest color bin. 
The large range of total IR properties in this bin suggests 
a mixture of galaxies with large amounts of dust emis- 
sion (LIRGs up to ULIRGs) and objects devoid of it. 
We note that, if we were to apply a different translation 
from MIR to total IR luminosity than simply averaging 
over the conversion factors derived from all reasonable 
templates, the observed trends remain intact. This is 
e.g. the case when an SED template corresponding to a 
larger heating intensity of the interstellar radiation fleld 
is used for objects with a larger rest-frame IR luminosity 
density, as done by Papovich et al. (2006; 2007). 

Since we divide our iiTs-band selected sample in bins 
containing equal numbers of objects, it is immediately 
clear that not only do red galaxies have on average the 
largest total IR luminosities, it is also true that they 
form the dominant contribution to the overall total IR 
luminosity emitted by ii's-selected galaxies at 1.5 < z < 
2.5. We show that the dependence of both color and Lir 
on stellar mass is an important driver of, but cannot fully 
explain the observed trends of Lir with color. 
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